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In this work we consider whether nonsymmorphic symmetries such as a glide plane can protect the 
existence of topological crystalline insulators and superconductors in three dimensions. In analogy 
to time-reversal symmetric insulators, we show that the presence of a glide gives rise to a quantized 
magnetoelectric polarizability, which we compute explicitly through the Chern-Simons 3-form of the 
bulk wave functions for a glide symmetric model. Our approach provides a measurable property for 
this insulator and naturally explains the connection with mirror symmetry protected insulators and 
the recently proposed TLi index for this phase. More generally, we prove that the magnetoelectric 
polarizability becomes quantized with any orientation-reversing space group symmetry. We also 
construct analogous examples of glide protected topological crystalline superconductors in classes D 
and C and discuss how bulk invariants are related to quantized surface thermal-Hall and spin-Hall 
responses. 


I. INTRODUCTION 

The last decade has seen a major breakthrough in the 
search of novel phases of matter with the discovery of 
topological insulators and superconductor^iM^. The orig¬ 
inal predictions of these systems have already led to many 
experimental realizations, in a very fruitful endeavor that 
continues today. The key insight underlying this dis¬ 
covery is that the presence of a symmetry, in this case 
time-reversal symmetry (T), allows to define a new bulk 
topological invariant of the Bloch wave functions in the 
Brillouin Zone (BZ). In a gapped fermion system, this 
invariant cannot change unless the gap closes, dehning 
a robust phase and protecting the existence of gapless 
boundary states. It was soon realized that other global 
symmetries in the Altland-Zirnbauer (AZ) classed, such 
as charge conjugation (C) and chiral symmetry, als o gi ve 
rise to new phases, leading to the periodic tabled of 
topological insulators and superconductors. 

The classification based on global symmetries then lead 
to the natural question of whether lattice symmetries 
can give rise to new topological phases of matter. For 
example, it was realized early on that in the presence 
of lattice translations, one may define extra topologi¬ 
cal invariants in lower dimensional slices of the Brillouin 
Zon^, which lead to the concept of weak topological 
insulator#*^] Point group symmetries such as rotations 
or reflections were also used to define new topological in¬ 
variants and phases of mat ter, w hich were termed topo¬ 
logical crystalline insulatorJi^^ (TCI) and superconduc- 
tors (TCSC). Recent efforts in the fieldi^fl^iT have been 
devoted to classifing these phases of matter protected 
by lattice symmetries in addition to global ones. Most 
of these previous works focused on symmorphic space 
groups for simplicity, i.e. groups where the full group is 
a semidirect product of the translation part and the point 
group. However, in view of the strong constraints that 
non-symmorphic symmetries place on the Bloch wave 
functions, one may expect that these symmetries can lead 


to richer structures, an idea that has drawn a lot of at¬ 
tention recently2lH2ll xhe question we consider in this 
paper is whether non-symmorphic symmetries, in partic¬ 
ular glide reflections, can define a new class of topological 
insulators or superconductors. We will focus on the three 
dimensional case without time-reversal symmetry, where 
it has been predicted that a new TCI protected by glide 
symmetry indeed exist^^. 

While in this work we will present explicit computa¬ 
tions of microscopic topological invariants, our main con¬ 
clusions can also be understood in a simple way by con¬ 
sidering topological bulk responses. It is well known that 
a three dimensional topological insulator can be char¬ 
acterized by a quantized bulk electromagnetic response 
term of the typ^^ 

which is known as the magnetoelectric response (because 
F A F cx E • B) or the “axion” Lagrangian. The second 
equality is expressed in coordinate-free notation, which 
we will use from now on (see Appendix]^. The magneto¬ 
electric coupling 0 is defined modulo 27r, and the presence 
of time-reversal symmetry requires that 0 = —0. This im¬ 
plies that 0 = 0 or TT, and the second case corresponds 
to a strong topological insulator. 6 can be computed mi¬ 
croscopically from the Chern-Simons 3-form of the Berry 
connection, establishing a direct correspondence with the 
^2 index. A physical consequence of 0 = tt is the pres¬ 
ence of an odd number of massless Dirac fermions on the 
surface which are protected by time-reversal symmetry. 

The same line of reasonin^^ implies that the magneto¬ 
electric coupling is quantized in the presence of any sym- 
metry that sends 0 —)■ —0, such as mirror reflectiorP^ISni 
which reverses one spatial coordinate. In the presence 
of a surface that respects this symmetry, one must also 
have an odd number of Dirac cones. A three dimensional 
topological insulator can therefore be protected by either 
time-reversal or mirror symmetry. The magnetoelectric 
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coupling in the second case can be computed microscop¬ 
ically from mirror Chern number^^ at invariant planes. 

The main result of this work is that glide symme¬ 
try also gives rise to quantized magnetoelectric coupling, 
which exactly corresponds to the TL^ invariant previ¬ 
ously defined for glide protected topological crystalline 
insulatord^. This result can be simply rationalized by 
the fact the magnetoelectric coupling is a bulk response 
property which makes no reference to the lattice. Since 
a glide differs from a regular mirror only by a half trans¬ 
lation, from the perspective of the bulk response both 
symmetries guarantee the quantization of the Q term, 
giving rise to a topological insulator when 9 = tt. This 
result is explicitly proven in Appendix [E| where we show 
in general that 9 becomes quantized in the presence of 
any orientation-reversing space group symmetry. 

In this work we also illustrate this general result by an 
explicit calculation of 9 via the Chern-Simons form for a 
particular model of a glide symmetric TI, conhrming the 
presence of a single Dirac fermion in the surface spec¬ 
trum when 9 = TT. In the second part of this work, we 
explain how these ideas can be naturally extended to su¬ 
perconductors without time-reversal symmetry in classes 
D and C. In these classes there are analogs of the mag¬ 
netoelectric coupling for the thermal and spin response, 
respectively, and these can also be quantized in the pres¬ 
ence of a glide symmetry. We will also present explicit 
models for these classes, and show that microscopic com¬ 
putations of the bulk topological invariants are consistent 
with the surface spectra. 


II. MAGNETOELECTRIC COUPLING AND Z2 
INVARIANT IN MIRROR AND GLIDE 
SYMMETRIC INSULATORS 

The quantized magnetoelectric response in Eq. has 
been long known to be th e dis tinguishing feature of 
strong topological insulator^^HM] time-reversal in¬ 

variance. The quantized coefficient in the action is known 
as the magnetoelectric polarization 9, and is given by the 
formula 

9 = — [ Tv(A^dA--A^A^J^. (2) 
47r Jbz V 3 y 

where = i (u”| dirt"*) for n and m conduction bands, 
and the integral is only gauge invariant modulo 27r, con¬ 
sistent with the ambiguity in Eq. While 9 is computed 
in this manner for translationally invariant systems, the 
corresponding topological response is robust against dis¬ 
order that preserves the symmetry on average and can 
be defined in the presence of interactions. The quanti¬ 
zation of 9 can also be prote cted by spatial symmetries 
like inversion symmetrySHMl improper rotationJi^, 

but cases where it is not possible to find a surface that 
preserves the symmetry lack protected gapless surface 
states. 


The existence of a quantized magnetoelectric response 
in a bulk material has an important implication for the 
response of surface states. In the presence of a small per¬ 
turbation that breaks time-reversal symmetry, the sur¬ 
face Dirac fermion becomes gapped, giving rise to a half¬ 
integer quantized Hall conductance^^. This cannot hap¬ 
pen in a pure two-dimensional system without topolog¬ 
ical order and reflects the topological nature of the 3d 
bulk. This behavior can be understood by considering 
the surface as an interface between the bulk and the vac¬ 
uum where 9 changes from tt to 0. The surface can thus 
be modeled by a spatially dependent 9 

such that 9{x < Q) = TT and 0(a; > 0) = 0 for an interface 
between the topological and trivial regions. Since F A F 
is a total derivative, one may integrate by parts to find 
the effective 2-1-1-D surface action of the Abelian Chern- 
Simons form 

S = M-f- [ A AdA (4) 

dTTh J 

which implies that the effective Hall conductance of the 
surface, given by A9 = tt in units of e^/(27r/i), is a half¬ 
integer value. This half-quantized topological response 
can serve as an additional feature to distinguish a topo¬ 
logical phase, and as we will see, can be generalized to 
other types of responses. 

To determine the coefficient 9 and identify a topolog¬ 
ical phase, one needs to explicitly evaluate Eq. The 
coefficient will be quantized in the presence of any sym¬ 
metry that takes 9 -a —9, but the computation may sim¬ 
plify in different ways for different symmetries. For ex¬ 
ample, in the presence of both time-reversal and inversion 
symmetry, Eq. can be related to the eigenvalues of the 
inversion operator at time-reversal invariant momentcP^. 
In the presence of a mirror symmetry, the computation 
can be related to the mirror Chern numbers at mirror 
invariant planes. In Appendix]^ we prove the quantiza¬ 
tion of 9 in the presence of a generic orientation-reversing 
space group symmetry with no reference to any of these 
simplifying circumstances. 

The purpose of this section is to demonstrate the ro¬ 
bust quantization of 9 imposed by mirror and glide sym¬ 
metries using microscopic tight-binding models. Starting 
with a brief review of crystalline insulators in class A with 
mirror symmetry, we show that the quantized magneto¬ 
electric response can be obtained as the parity of the the 
integer-valued (S) topological invariant protected by the 
mirror. Howev er, when the mirror symmetry is replaced 
by a glid^22Illl^ g^ly a ^2 invariant survives, which cor¬ 
responds to the quantized 0 = 0, tt. 
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A. Mirror symmetry 


B. Glide symmetry 


In three dimensional insulators with mirror symme¬ 
try, topological invariants can be defined by considering 
the mirror invariant planes in the BZ, where bands have 
a definite mirror parity. Total and mirror Chern num- 
bers can be thus defined for these invariant planeP^lMl 
The total Chern numbers for cuts perpendicular to the 
mirror plane vanish by symmetry, and for simplicity we 
assume they are zero in the mirror-invariant planes as 
welP^. In a system with mirror symmetry reflecting the 
z axis, the bulk BZ has two mirror invariant planed 
at /cz = 0 and tt (Fig. (a)). This allows us to label 
bands in these two planes by their mirror eigenvalue 
(F = 0 and 1 for spinless and spinful fermions respec¬ 
tively), as no terms mixing the two sectors are allowed 
by symmetry. Chern numbers for the even and odd 
occupied bands are separately well defined for = 0 
and TT. The mirror Chern number for a mirror-invariant 
plane {k^ = 0, tt) is defined as the difference between the 
two sectors ^ — C^T ). Consider, for example 

the case with nonzero Chern numbers for the even and 
odd sectors in the kz = 0 plane {Cq = ±1) and vanishing 
Chern number for both sectors at kz = tt (C^ = 0), now 
= 1 and = 0. 

A minimal Hamiltonian implementing this phase can 
be obtained starting from t he 4-b and model of a 3d TI 
with time-reversal symmetr^ F^FS I 7 ~ = iayJC: 


= tx sin kxTy + ty sin kyazTx + tz sin kzCTyTx -b mkT^ 

( 5 ) 

with TTT-k = m — cos ky. This model is in the strong 
TI phase for I < \m\ < 3 with 9 = tt. Now we can re¬ 
move the time-reversal symmetry constraint and instead 
demand invariance under mirror M = i^az reflecting 
the z direction. The Hamiltonian ([^ formally has both 
symmetries, but the mirror allows different perturbations 
than time-reversal. It is easy to see that the two mirror 
sectors (with opposite az eigenvalues) now have oppo¬ 
site Chern numbers in the kz = 0 plane and vanishing 
Chern numbers at kz = tt . Similar to the weak indices 
in time-reversal invariant insulators, and sepa¬ 
rately rely on translational symmetry along z direction 
but the strong mirror Chern number -I- is 

robust against translational symmetry breaking, as long 
as mirror symmetry is preserved. We note that, with 
vanishing total Chern number {i.e. = 0), the 

quantized magnetoelectric coupling 9 is completely de¬ 
termined by the strong mirror Chern number, 9 = ttC^^ 
(mod 2tt). This is seen by counting surface modes: in 
the invariant planes = —C^ counts the number of 
chiral modes on the surface at kz — 0 oi tt propagating 
right (left) in mirror sector -|- (—). Each pair of counter- 
propagating modes forms a surface Dirac-cone, so if the 
total number is odd, the bulk has nontrivial 9. 


In this section we present an alternate picture to un¬ 
derstand the glide-protected TL’i invariant first proposed 
in Ref [2SI (see also Appendix]^, in terms of the quan¬ 
tized magnetoelectric polarization. Let the glide G = 
{Mz|ta;/2} reflect the z direction and translate along x 
by half of a unit cell (Fig. [^. = {(—I)^l|ta;} is a 

pure translation, with the sign depending on how a 27r 
rotation is represented. The eigenvalues of the glide op¬ 
erator are As we traverse the Brillouin Zone 

(BZ) in the x direction on a line that is pointwise invari¬ 
ant under the symmetry, the eigenvalues wind into each 
other. In a system respecting this symmetry if we follow 
a band with the eigenvalue, it is connected to 

band with at the zone edge, so the bound¬ 

ary condition for the Bloch wave functions is constrained 
to |uk_|_G ) I'^ic) where is the reciprocal lattice 
vector parallel to tj,. So bands in the presence of glide 
symmetry come in pairs that cannot be separated by a 
gap, as the crossings are protected by the fact that the 
pairs have different eigenvalues under the symmetry. The 
Chern number for each single band is ill defined as one 
band must evolve into the other on the zone boundary, 
only the total Chern number for the pair is a topological 
invariant. 

Again, for simplicity, we assume that the total Chern 
number for conduction/valence bands in any 2d cut par¬ 
allel to the mirror plane vanishes. Nonzero values for 
perpendicular cuts are forbidden by mirror symmetry. 
The minimal model realizing the nontrivial phas^^ is 
analogous to the mirror-symmetric case i.e. Eq. ([^ 




AG 


tx sin 



PxTx + ty sin kyTyF 


-I- tz sin kzPzTx + rai^Tz 


( 6 ) 


with the glide operator Gk = F'^^l'^px- Here Wk has 

the same form as in ([^. It is easy to see that the model 
is gapped for appropriate choice of parameters and band 
degeneracies can be removed almost everywhere in the 
BZ with symmetry allowed terms. Pauli matrices t and 
p act on orbital and sublattice degrees of freedom re¬ 
spectively and the Hamiltonian preserves glide symmetry 
G\iH\^ = with Mz {kx^ ky, kz^ = {kx^ky, kz)- 

We used the convention where operators and Bloch wave 
functions are not periodic in the BZ (see Appendix [d| , 
and the model can be regarded as either spinless {F = 0) 
or spinful {F = 1) with spin-polarized electrons, such 
that trivial bands with opposite z-component spin 
are pushed far over or below the Fermi level and can be 
omitted (note that Mz is diagonal in the basis). 

Regardless of the microscopic differences, in the macro¬ 
scopic translational-invariant electromagnetic response 
theory there should be no distinction between a mirror 
and a glide, and, as we formally prove in Appendix 
9 is quantized as 0 or tt just like with a mirror. We 
numerica lljEH verified that in the nontrivial phase of this 
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FIG. 1. (a) Bulk BZ of a glide or mirror symmetric crystal 
with the two invariant planes BZ in red (shaded) and invariant 
lines in blue (thick) if particle-hole symmetry is also present, 
(b) Surface BZ for a cut normal to y with invariant lines in 
red (thick), labels for high symmetry points and a sketch of 
an occupied band pair along these lines. On the invariant 
planes {k^ = 0 and kz = it) the glide eigenvalue (color code) 
distinguishes the two bands, (c) Example of a crystal with 
glide symmetry G that reflects the 2 direction and translates 
by half of the unit cell in the x direction, the glide operator 
exchanges the two sublattices (A and v)- This pattern is 
repeated in parallell planes shifted perpendicular to the plane 
of the drawing. 


FIG. 2. (a) Band structure in slab geometry with 20 unit cells 
in the y direction for glide protected class A topological insu¬ 
lator. We use with = 1, m = 2 and (j> = 0.4. The sur¬ 
face Dirac cone is at a generic momentum on the high symme¬ 
try line, the left and right moving branches are distinguished 
by their glide eigenvalues and cannot gap out. (b) Evolution 
of 0 from numerics while tuning across the transition from 
the topological to the trivial phase without breaking the glide 
symmetry. The error bars indicate two standard deviations 
(95% confidence interval) of the Monte Carlo estimates. 


a mirror or a glide. We now show how these consider¬ 
ations can also be applied to superconductors in three 
dimensions where T is broken. In this section we con¬ 
sider a superconductor with no other local symmetry but 
the particle-hole symmetry of the Bogoliubov-de Gennes 
Hamiltonian, which belongs to class D. We will show that 
the presence of an additional glide symmetry protects the 
existence of a topological crystalline superconductor with 
a single Majorana cone at the surface. In section [TV] we 
will consider the analogous problem in the presence of 
SU(2) spin rotation symmetry. 


model, 0 = TT is robust against symmetry preserving per¬ 
turbations (Fig.[^ for details see Appendix [A|. We would 
like to point out that starting from a mirror symmetric 
TCI (for example 0) one can double the unit cell and 
weakly break down separate mirror and half translation 
symmetries to a glide without closing the gap or chang¬ 
ing the value of 9. While the mirror Chern numbers are 
no longer well defined, the TLi invariant defined by 9 sur¬ 
vives. We emphasize that 0 is a macroscopic response 
quantized by macroscopic mirror symmetry, so it is ro¬ 
bust against symmetry preserving interactions and dis¬ 
orders which preserve the symmetry on averag^^, both 
in the case of glide and mirror symmetry. To summarize, 
we provided a physical understanding of the TL-i invariant 
introduced in Ref. [231 in terms of quantized magnetoelec¬ 
tric coupling. 


III. TOPOLOGICAL CRYSTALLINE 
SUPERCONDUCTORS IN CLASS D 

In the previous section we have shown how the mag¬ 
netoelectric coupling in insulators can be quantized in 
the presence of symmetries other than T, in particular 


A. Bulk invariant and surface thermal Hall 
conductance 

The reason why a glide can protect a topological super¬ 
conductor in three dimensions without time reversal sym¬ 
metry is that this phase is also characterized by 0 term 
that is quantized with any orientation-reversing symme¬ 
try. The reasoning is analogous to the one used for in¬ 
sulators. We first consider the case of a superconductor 
with time-reversal symmetry, in class DHL In three di¬ 
mensions, class Dill has an integer topological invariant 
1 / S S which counts the number of Majorana cones at the 
surface. In the same way as the insulator, in the presence 
of a weak perturbation that breaks time reversal symme¬ 
try, the surface becomes gapped and each Majorana cone 
contributes half of the minimal thermal Hall conductance 
of a 2d supercond uctor (that is half of the minimal value 

for a 2d insulator (pSEHl 

l^xy _ (tt/cb)^ V 

T “ 3h 4' 

Formally, a class D superconductor is the same as an 
insulator with an extra antiunitary particle-hole sym- 
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metry C that anticommutes with the Hamiltonian and 
squares to +1, because a particle-hole symmetric Bloch 
Hamiltonian for insulators has the same form as the 
Bogoliubov-de Gennes (BdG) Hamiltonian for supercon¬ 
ductors. The only important difference is that only half 
of the degrees of freedom in the BdG Hamiltonian are 
physical, since all negative energy states correspond to 
the annihilation operators of the positive energy excita¬ 
tions over the BdG ground state. This is why the surface 
Dirac cones in insulator case reduce to surface Majorana 
cones in the superconductor case. 

The fact that the insulator and superconductor prob¬ 
lems are formally the same allows us to use the BdG 
Hamiltonian the same way as the Bloch Hamiltonian to 
calculate 9 from the band structure of a glide-symmetric 
superconductor, which must be quantized to 0, tt by the 
same reason as in the insulator case. When 0 = tt, this 
implies an odd number of surface Majorana cones, and 
a half-integer thermal Hall conductance when the glide 
symmetry is broken. This is not allowed in a purely two- 
dimensional gapped superconductor with no ground state 
degeneracy, where K^y/T is always quantized in integer 
multiples of {-kUbY/(& h). 


(a) (b) 



FIG. 3. Band structures in slab geometry for glide protected 
class D topological superconductor. We use with ty = 

— 1. (a) with m = 2.5 the surface Majorna cone is at M, while 
(b) with m = 0.5 at X. The left and right moving branches 
are distinguished by their glide eigenvalues on the invariant 
lines and cannot gap out. Inset: Evolution of 9 while tuning 
across the transition from the topological (m < 3) to the 
trivial (m > 3) phase from numerics. Note that at m = 3 the 
gap closes and 9 takes on an intermediate value not allowed in 
a gapped system. The error bars for the Monte Carlo results 
are smaller than the symbols. 


C. Lower dimensional topological invariants 


B. Microscopic model with glide plane in 3d 


The explicit model for a glide symmetric superconduc¬ 
tor in class D is very similar to the insulator in Eq. Q, 
but with an extra particle-hole symmetry = 1. Again 
we consider a glide plane reflecting the z direction and 
translating along x, G = {Mz\tx/2}. We represent 
charge conjugation as C = t^JC where JC is the complex 
conjugation operator and the r act on the particle-hole 
degree of freedom. As charge conjugation acts locally, it 
has to commute with the glide as G_kC = CGk- The 
Hamiltonian with (f> = 0 possesses these symmetries 
for spinless fermions, however, such systems do not nat¬ 
urally appear as superconductors and require fine tuning 
as insulators. For the rest of this section we assume the 
physical case of spinful fermions, for which the appropri¬ 
ate choice of the glide representation is Gk = ie^^^^^pxTz 
(see Appendix]^ for details about symmetry representa¬ 
tions in BdG systems). A simple Hamiltonian respecting 
these symmetries is: 


k 

Hk — ix cos Py^x fy sin kyPzTx tz 


Cl n Zr* 


For appropriate choice of parameters (see Fig. the 
model is in its topological phase with gapless excitations 
on symmetry preserving surfaces and numerical evalua¬ 
tion confirms that 6 = tt (Appendix[A|). Fig.|^(a) and (b) 
shows surface spectra with single Majorana cones pinned 
by particle-hole symmetry to different high symmetry 
points. As will be explained later, in fact the Majorana 
cone can only appear at X or M, but not T or Z in 
the surface BZ (shown Fig.[^(b)). This is an important 
difference compared to the insulator case. 


The structure of the surface modes in the partic¬ 
ular case of class D can be further understood from 
the presence of lower dimensional topological invariants 
(known as weak invariants or indices) associated to glide- 
invariant lines and planes in the Brillouin Zone. While 
our model has been chosen such that all Id and 2d in¬ 
variants associated to particle-hole symmetry are trivial, 
the presence of a mirror or glide enables new lower di¬ 
mensional invariants. 

First we review the case of mirror symmetrjff^E^. The 
classification for TGSGs in class D depends on the square 
of mirror operator i.e. = ±1. In particular, no 

nontrivial TGSGs exist with M_, this means, while 9 is 
quantized, only 0 = 0 is allowed by symmetry. On the 
other hand, with the choice of M_|_, we find a rich struc¬ 
ture with both Id and 2d weak invariants. Similarly to 
class A, there is an integer-valued 2d index (mirror Ghern 
number) in mirror-invariant planes. Besides, there are Id 
mirror S 2 invariants along the high symmetry bulk lines 
in y direction (blue lines in Fig. (a)). These mirror 
S 2 invariants guarantee the presence of a pair of sur¬ 
face zero modes of opposite parity at the corresponding 
surface high-symmetry momentum, which split for any 
other momenta generating a Majorana cone. Therefore, 
these Id mirror indices determine the parity of the num¬ 
ber of surface Majorana cones on the surface projections 
of the mirror invariant lines (Fig.[^(b)). This leads to a 
TL\ index in the case of M+ mirror symmetry, as long as 
translation symmetry is preserved. 

The case of a glide can be understood as a combina¬ 
tion of the above two cases as we argue below. Mirror 
Ghern numbers cannot be dehned with a glide for the 
same reason that applies in class A. However, the smaller 
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symmetry group with glide symmetry poses weaker con¬ 
straints on the band structure compared to mirror M_, 
allowing 0 = tt as illustrated by our model. Id indices 
can still be defined, but the square of the glide operator is 
different for high symmetry points since it changes with 
kx, unlike for a mirror. For the lines a,t kx = 0 we have 
Gq = —1, but for kx = TT, = 1. This diffe rence is key 
because for Gq = —1, no mirror index exist^i^E^, while 
for G^ = 1 there is a S 2 index at each high symmetry 
point. Therefore, surface Majorana fermion can only be 
found at high-symmetry points with kx = tt, i.e. the X 
or M points. Two cases with a single Majorana cone at 
X OT M are realized in our proposed model for m = 0.5 
and m = 2.5, as seen in Fig. in both cases 9 = tt 
as it is determined by the parity of the total number of 
Major anas. 

D. Surface Dirac model 

An alternative approach to demonstrate the protection 
of a topological phase by a symmetry is to consider how 
symmetries are implemented in a generic surface theory. 
For example, for a regular topological insulator, the pres¬ 
ence of time-reversal symmetry protects a single Dirac 
cone to be gapless. If a single Dirac cone is found at the 
surface, it cannot be removed until the symmetry is bro¬ 
ken or the bulk gap closes. Two surface Dirac cones can 
however be gapped without breaking the symmetry. 

We consider how the glide symmetry is implemented in 
a generic Dirac Hamiltonian at a high-symmetry point, 
offering an alternative explanation of the results in the 
previous section. As mentioned earlier, as a result of 
particle-hole symmetry a single Majorana cone can only 
appear at a high symmetry surface momentum, we will 
discuss the case of more Majoranas later. In a 2-band 
model Go that squares to —1 can be chosen iai for 
i = x,y,z, but since ax and behave the same un¬ 
der complex conjugation we only need to consider ax or 
ay. Charge conjugation is represented as C = UqJC with 
real unitary Uc ■ A 2-band gapless Dirac Hamiltonian has 
the form 

H = kxTx + kzT^ ( 9 ) 

where the F’s are hermitian, anticommuting, with ±1 
eigenvalues and [Fa,, Go] = {r 2 ,Go} = [ri,C] = 0. If we 
choose Go = iay, [Uc,Go] = 0, so Uc = 1, Fa, cannot be 
chosen to satisfy all commutation and anticommutation 
relations. Similarly if Gq = iax, {Uc,Go} = 0, Uc = a^ 
and again no Fa, is allowed. Therefore it is impossible to 
write a Dirac Hamiltonian with particle-hole symmetry 
and Gq = —1. This shows that a single surface Majorana 
cone is forbidden at the F and Z points of the surface BZ 
(see Fig. [^(b)) in the presence of a glide. 

On the other hand with G^ = -|-1 we can choose 
Gt, = ax, Uc = Tx = ax and = a^- Now a sin¬ 
gle Majorana cone is allowed, but it cannot be gapped 
out, for that we would need a mass term toFo such that 


{ro,ri} = {FqjC} = [ro,G,r] = o. One can check that a 
mass term is not allowed for any valid choice of a 2 x 2 
representation, single Majorana cones are allowed and 
protected at X and M. 

Finally, we may also consider a system with a pair 
of cones at opposite surface momenta. Similarly to class 
Aim, 

a pair of surface Majorana cones with different glide 
eigenvalues at one high symmetry point are locally pro¬ 
tected, but can symmetrically move around the BZ and 
gap each other out at another point where their eigen¬ 
values are the same. This shows that only the number of 
cones modulo 2 at each of X and M is stable against 
symmetry-preserving perturbations. The classification 
with full translation invariance is thus , while allowing 
terms doubling the unit cell in the 2 direction reduces 
the classification to TLi counting the parity of the total 
number of surface Majorana cones. Such a S 2 index is 
given by the Chern-Simons 3-form in Eqn. (©■ 

IV. CLASS C SUPERCONDUCTOR WITH 
GLIDE PLANE IN 3D 

In this last section we consider how a singlet super¬ 
conductor with SU(2) spin rotational symmetry in three 
dimensions may also have a topological phase. This type 
of superconductor belongs to class C. After appropriate 
rearrangement of the degrees of freedom (see Appendix 
0 the Hamiltonian can be block-diagonalized in the spin¬ 
s'"^ basis where the two blocks are unitarily related with 
identical spectra and topological properties. In the re¬ 
duced problem charge conjugation is combined with a 
spin rotation acting as C = r^/C, with = —1, which is 
the main difference compared to class D. 

To understand the emergence of the topological phase 
protected by a glide in class C, it is instructive to first 
consider the more familiar case of an SU(2) invariant 
superconductor with time-reversal symmetry, which be¬ 
longs to class Cl, and the anomalous response of the sur¬ 
face after breaking time reversal symmetry. After this, 
we argue how the 9 term is also quantized in a class C 
superconductor with a glide, present an explicit model 
for this, and also argue how the protection of the surface 
states can be seen directly from the surface theory. 

A. Bulk invariant and SU(2) axion term 

In class Cl, topological superconductors are charac¬ 
terized by a topological invariant v, which counts the 
number of pairs of surface Majorana cones. The SU(2) 
spin rotation symmetry allows the definition of spin Hall 
conductances]. Once the surface is gapped by breaking 
time-reversal symmetry (but not SU(2) spin rotational 
symmetry), the surface spin quantum Hall conductance 
is given bj®* 

_ mf V 


h 2 


( 10 ) 
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where v is an even integer. Note that when v = 2 the 
above a^y is only half of that of a d + id singlet super¬ 
conductor in 2cP. Therefore this anomalous half-integer 
surface spin quantum Hall conductance serves as a probe 
to characterize the nontrivial 3d phase. It should be 
noted that the thermal Hall conductance Kxy is well de¬ 
fined, but it is not sufficient to characterize the topolog¬ 
ical phase. 

This half-integer response can be related to an analog 
SU(2) bulk 9 term, in an analogous way to a 3d topolog¬ 
ical insulator. To see this, we first consider the effective 
SU(2) continuum gauge theory that describes the spin 
quantum Hall superconductor in 2-1-1 dimensions, which 
captures the response of the SU(2) spin rotation invari¬ 
ant system coupled to a SU(2) gauge field. While this 
gauge field is fictitious, this treatment is useful to derive 
the response to an external Zeeman fielcP^. This system 
is described by the effective actioiP^ 

S ^ f AdA^ ( 11 ) 

dn h J 

where A^ is the z component (in spin space) of the SU(2) 
gauge field, which we identify as the z component of an 
external Zeemann field A^ = . C = f B € 2Z is the 

Chern-number of the negative energy bands in one spin 
sector, which is an even integer in class C. This is analo¬ 
gous to the quantum Hall effect, but as we are interested 
in spin currents, the electric charge e is replaced by /i/2 
in the coupling (we take the g-factor g = 1). The spin 
current is 


J* = SS/SAf = (12) 

so the spin Hall conductance is a^y = C. To get 
a proper SU(2) gauge theory, we promote H to a non- 
abelian SU(2) gauge field hy A = the action com¬ 

patible with the previous one is given by the nonabelian 
Chern-Simons 3-form 


1 (hl2f G 
Att h 2 


Tr ( H A dH - y H A H A H ) , (13) 


where an extra 1/2 appears to compensate for the trace. 
As we see, compared to the U(l) gauge theory, the co¬ 
efficient is only half of the Chern number for one spin 
sector. 

In analogy to Eq. Q, an SU(2) axion action can be 
defined in 3-1-1 dimensions in the following way: 


1 (fi/ 2)2 6 » 

Att h 2 


TtF A F 


(14) 


where F is the nonabelian field strength tensor and 9 is 
the Berry Chern-Simons 3-form (Eq. (§) for the occu¬ 
pied bands in one spin species. The ambiguity in 9 is 
now 47r because of the extra factor of 1/2. One can also 
check that a spatial domain wall of 47r in theta gives rise 
to a surface with the minimum allowed value of the spin 


Hall conductance in a 2d system, corresponding to C = 2 
in one spin sector. As we argued in class A, if the sys¬ 
tem has a symmetry on average that flips an odd number 
of spacetime dimensions, 9 is quantized to 0 or 27r mod¬ 
ulo 47r, leading to a 2^2 classification. The reason why 
the ambiguity in 9 calculated from the band structure 
changes to 47r is quite subtl^^. It comes from the gauge 
fixing requirement that the band structure can be con¬ 
tinuously deformed to the trivial band structure in the 
trivial gauge without breaking C. We refer the interested 
reader to Appendix for details. 


B. Microscopic model with glide plane 


Following the same logic as for classes A and D, we 
now consider how the presence of a glide symmetry can 
protect an SU(2) invariant topological superconductor in 
the absence of time-reversal symmetry, i.e. in class C. 
As in the previous classes, the presence of an orientation 
reversing symmetry is sufficient to guarantee the quanti¬ 
zation of the 9 term (now to either 0 or 27r), which is the 
topological invariant that characterizes the phase. This 
phase has formally the same properties as a Cl topolog¬ 
ical superconductor with v = 2, namely a pair of Majo- 
rana surface cones and a half-quantized spin Hall conduc¬ 
tance upon breaking the glide symmetry on the surface. 

To show this, we now consider a microscopic model in 
class C with a glide symmetry, demonstrating the pres¬ 
ence of protected surface modes and computing the value 
of 9 explicitly. In the original full Hilbert space the natu¬ 
ral representation of the glide is Gk = Px^zC^z-, but 

because of the full SU(2) spin symmetry we can cancel 
the spin rotation part by attaching —iTzUz to our defi¬ 
nition, so we may use Gk = F'^^/'^px- For this operator 
G^ = -|-e®^“= 1, showing the “spinless” nature of the prob¬ 
lem. A Hamiltonian for one spin component with these 
symmetries can be constructed as 

k 

=Axy sin ky sin -^PxTx + Axz sin kx sin kzPzTx+ 

+ Ao(cos kx — cos ky + a)Ty + m^Tz- (15) 

For appropriate choice of parameters (Fig. this 

Hamiltonian realizes a nontrivial topological phase with 
an odd number of pairs of surface Majorana cones (not 
counting the spin degeneracy). Class C is similar to class 
A as it has a 2TL index in 2d without symmetry, corre¬ 
sponding to an even total Chern number in xy cuts which 
vanishes in our model and we restrict our further discus¬ 
sion to this case. Again, no integer-valued mirror Chern 
number can be defined in glide-invariant planes, but a 
2^2 index still survives due to the bulk quantization of 
0, which we have computed explicitly (Fig. |^inset). 








FIG. 4. Band structure in slab geometry for glide pro¬ 
tected class C topological superconductor. We use with 
Axy = A^cz = Ao = 1, a = 0.5 and m = 1.5. Note the 
symmetric pair of surface Majorana cones at a generic mo¬ 
mentum, the left and right moving branches are protected 
from gapping out by different glide eigenvalues. This is the 
spectrum for one spin sector, for the full system there is an 
additional spin degeneracy. Inset: Evolution of 9 while tuning 
across the transition from the topological to the trivial phase 
from numerics. Note that at m = 2.5 the gap closes and 9 
takes on an intermediate value not allowed in a gapped sys¬ 
tem. The error bars for the Monte Carlo results are smaller 
than the symbols. 


C. Surface Dirac model 

The fact that an odd number of pairs of Majorana 
cones is protected in class C with a glide can be shown 
explicitly by considering the surface Hamiltonian of a 
pair of Majorana cones with = —1. It is instructive to 
show first the known case of class CL The Hamiltonian in 
the vicinity of a high symmetry surface momentum (Fig. 
0(b)) is 

H = axK + cTzkz: (16) 

where the Dirac matrices are 4x4 and spanned by cji, r^, 
the identity is implicit. The particle-hole operator is C = 
iTylC and the time-reversal operator is T = UyTyK, with 
K complex conjugation. There are four possible mass 
terms for this Hamiltonian, ayTx, <yyTy, (TyTz, ay. The first 
three masses are forbidden by C, and the last one, ay, is 
forbidden by T. Therefore, a single pair of Majorana 
cones cannot be gapped out at the surface in class CL 
In class C, when time reversal symmetry is broken, 
both a mirror or a glide can still protect the presence of 
a single pair Majorana cones, because as we now show, 
both of these symmetries forbid the ay mass as well. In 
the presence of a mirror, we consider a reflection z —)■ 
—z, with operator M± that satishes = ±1. In the 


presence of a glide, we would have but since 

the Majorana cones must be at opposite momenta due to 
C, their annihilation can only occur at = 0,tt and we 
only need to consider these two cases with Gq = -1-1 and 
Gj = —1, same as the situation of regular mirrors. We 
now discuss the two cases of M±, which apply to both 
mirror and glide. 

In the M+ case, for the two eigenstates of ax with 
positive eigenvalues, the mirror eigenvalues are the same 
and equal to 1, while the other two are equal to -1, as 
would happen if M were the glide operator at fc = 0. The 
operator doing this is simply M+ = ax- In the M_ case, 
for the two eigenstates of ax with positive eigenvalues, the 
mirror eigenvalues are , as would happen if M were 
the glide operator aX k = tt. This is achieved with M_ = 
iaxTz- Note both M± anticommute with because M 
reverses kz- Also note both M± commute with C as we 
want. 

Once we have the operators for M±, it is easy to see 
that both operators forbid the mass ay. Since the other 
three masses are already forbidden by C, a pair of Majo¬ 
rana cones is protected in the presence of = —1 and 
the additional mirror/glide symmetry M±. 

V. CONCLUSIONS 

In this work we classified 3d topological insulators 
and superconductors protected by non-symmorphic glide 
symmetry in classes A, D and C and presented lattice 
models for these phases. Our results, however, are more 
general. As our arguments only rely on symmetries of the 
effective long-wavelength response theory, the S 2 classi¬ 
fication is also robust if mirror symmetry, or any sym¬ 
metry reversing an odd number of spatial coordinates, 
is preserved on average. The cases with glide symme¬ 
try illustrate that these macroscopic considerations are 
insensitive to the fractional translation that accompa¬ 
nies the mirror operation, and identify the most robust 
topological invariants that are also defined (among oth¬ 
ers) with simple mirror symmetry. Glide symmetry is 
present in over a hundred of the 230 crystallographic 
space groups, and all but the 65 chiral groups contain 
orientation-reversing operations, so our results should be 
widely applicable to experimentally and numerically de¬ 
tect topological crystalline insulators and superconduc¬ 
tors without time-reversal symmetry. 
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Appendix A: Evaluation of Chern-Simons 3-form 
and second Chern-form 

In this section we review the calculation of the magne¬ 
toelectric coupling. All our formalism is in k-space, d de¬ 
notes the exterior derivative and A the wedge product!^, 
these are antisymmetrized in the spatial indices (/i, v, 
...). Operators either act on the occupied band space 
and trace is taken over occupied band indices (n, m, 
...) or the full Hilbert space of the unit cell and we 
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supress most indices for the sake of a compact notation. 
The magnetoelectric coupling is defined in terms of the 
Chern-Simons 3-form: 

0 = — [ T^:(A^AA--A^A^J\^. (Al) 
47r Jbz V 3 y 

The main difficulty about evaluating this expression is 
that one has to find a (patchwise) smooth gauge in the oc¬ 
cupied band space, which is a complicated task in numer¬ 
ical studies. To circumvent this, we instead find a gapped 
deformation to a trivial state with constant Hamiltonian 
using a tuning parameter ^ 4 , such that 9{ki = 0) = 0, 
and calcula te the change in 0 by the 4 dimensional second 
Chern forrEpnEH that is locally gauge invariant: 

0{7t)- 9{0) = ]- [ r Tt{AAA) (A2) 

47’' Jbz Jki=o 

where J- = AA — iA A A, the nonabelian Berry curva¬ 
ture in the 4 dimensional space spanned by = (k, k^). 
(Note that in the convention we use J- A A = 
We realize that IF can be written in 
a gauge invariant form as J" = iV{AV) A {AV)V using the 
projector onto occupied bands Vy = J2neocc. Wk) (^k L 
which is also gauge invariant. In this formulation is an 
operator acting on the full Hilbert-space of the unit cell, 
but only nonzero on occupied bands, so we can extend 
the trace to the full Hilbert-space without changing the 
result. The usual components can be obtained as ma¬ 
trix elements between occupied states in a given basis, 
^k'^ = {u^\Fk \ u'^)- So we arrive at the locally gauge- 
invariant expression 


6»(7r) -6»(0) = 


(A3) 




Tr {V{AV) A [AV)V A [AT) A {AV)) 


that we numerically evaluate using adaptive Monte Carlo 
integratioiP^. In order to get physical result, one must 
use the Bloch formalism where the orbital positions are 
taken into account ( Appendix|^and|E ), the interpolating 
Hamiltonians in classes A and D can be written as 

=\{l-cosk,)H^^+ 

-\- sin /c 4 sm{kx/‘J)py + ^ (1 + cos k 4 )Tz (A4) 

Hk° =^(1 - coski)H^^+ 

-\- sin ki sm{kx/‘J)pxTx + ^ (1 + cos k4)Tz (A5) 

with and given in ^ and (j^ respectively. We 
find that 9 is quantized to 0 or tt (mod 27r) to high ac¬ 
curacy, a result that is robust against symmetry allowed 
perturbations of the final Hamiltonian and deformations 
of the interpolation as long as the bandgap does not close 
(Fig.[|(a)). 



k, k, 


FIG. 5. (a) Evolution of 6 during the gapped deformation 
from the trivial to the topological phase while breaking 
glide symmetry in class A. We use parameters ty, — 1, 0 = 0.4 
and m = 2 (circles) or m = 2.5 (squares) and we also show 
a deformation with m = 2 but with the symmetry allowed 
perturbation ^Ty with fi = 0.5 added to (triangles), (b) 
Same plot for class C using with /A^y = = Aq = 1, 

a = 0.5 and m = 1.5 (circles) or m = 2 (squares) and m = 1.5 
with symmetry allowed perturbation Atx with /3 = 0.5 added 
to . The error bars are smaller than the symbols. 


In class C we have to reconcile the change of am¬ 
biguity in 9 from 27r to 47 r from the band structure 
point of view. A natural gauge choice for a class C 
system is the requirement that = iTy that 

relates unoccupied states at k with occupied states at 
—k. Gauge transformations preserving this condition 
satisfy Tyll^ = {U't^*Ty where [/“/“ act on the occu¬ 
pied/unoccupied bands. This constraint is not sufficient 
to remove the 27r gauge ambiguity in 9 coming from the 
winding number of U^, as any is allowed as long as it 
is accompanied by the appropriate U^. 

So how can we define a proper bulk invariant? The 
idea is to prove that for a cyclic gapped deformation of 
the band structure (with k^ as the tuning parameter) the 
second Chern form 


1 

47r 


TyAAA 


(A6) 


is quantized to multiples of 47 r as long as particle-hole 
symmetry in any 3d cut is preserved (as opposed to multi¬ 
ples of 27r without symmetry). This is proved in Ref. HU 
for our case H = 3 is the dimension of fc-space and <5=1 
the dimension of real space, for us this is fc 4 , the tuning 
parameter that is not affected by particle-hole symmetry. 
The result is, for D — S = 2 in class C the classification 
is 2Z, proving our conjecture. 

Now, if we find a gapped particle-hole symmetric de¬ 
formation from the trivial band structure to the glide 
symmetric one, we can calculate the difference A9 = 
9{'k) — 0(0) in the 3d Chern-Simons forms between the 
initial and final states using the locally gauge invariant 
expression of the second Chern form in terms of the band 
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projector. 27r will be different from zero, as the ambiguity 
introduced by different deformations is 47r. Such defor¬ 
mation to the trivial state always exists, as there are no 
nontrivial phases in 3d class C without any other sym¬ 
metry. Note that the formula in terms of T requires a 
continuous gauge choice throughout the deformation and 
0(0) is only zero in the trivial gauge. We can view this as 
a nontrivial gauge fixing condition demanding that the 
band structure is continuously deformable to the trivial 
state in the trivial gauge, in this gauge the Chern-Simons 
3-form in the final state gives the same result. However, 
we know of no method of checking whether this condition 
is satisfied other than explicitelly constructing a deforma¬ 
tion. On the other hand, the formula with the projector 
is completely gauge-invariant, insensitive to discontinu¬ 
ous large gauge changes that would add extra 27r’s to the 
formula with so we do not have to worry about the 
band structure satisfying any gauge condition using this 
method of computation. 


For the specific model in equation (15) we use the de¬ 
formation preserving C = TylC : 


H, 


CG 


= x(l — COSfc4)iJ^'^-|- 


1 


-I- siTik4sin{kj;/2)pyTj; -|- -(1 -I- cos/c4)r2. (A7) 


Our numerical results give 0 = 0 and 27r (mod 47r) 
to a high accuracy for the trivial and topological phases 
respectively (Fig. Kb)). 


Appendix B: Relation to earlier definition of K 2 
index in class A 

Here we review the alternate definition of the TL-i index 
in class A introduced in Ref. [231 There are two planes 
in the 3d BZ that are invariant under a mirror or glide, 
at fcz = 0 and tt, in the rest of the BZ glide does not act 
locally, bands do not have a well defined eigenvalue and 
symmetry allowed perturbations generically destroy most 
degeneracies. A 2 x 2 Hamiltonian without any symme¬ 
try constraints only has pointlike degeneracies between 
occupied bands in 3d (Weyl nodes), which allows us to 
choose a surface in the BZ that connects the k^ = 0 and tt 
planes such that the all the bands are nondegenerate. We 
will choose a constant k^ surface for simplicity and show 
later that the result is insensitive to this choice or the 
assumption that there are only pointlike degeneracies. 

Now we can choose a surface in the bulk BZ that in¬ 
cludes the two invariant planes plus a surface connecting 
them traversing half the BZ, say at = 0 (see Figj^(a)). 
We can follow a band around this surface (Figj^(b)), this 
is possible because of the well defined glide eigenvaue in 
the invariant planes and the generic lack of degenera¬ 
cies in the connecting surfaces. We can compute the 
“glide Chern number” corresponding to one such band 
on the “bent BZ’^^ by counting the winding number of 




FIG. 6. Illustration of our definition of the glide Chern num¬ 
ber. (a) Bulk BZ of a glide symmetric crystal with a possible 
choice of the “bent” BZ in red (shaded). The glide Chern 
number is calculated by integrating the Berry flux following 
a single band around this torus, (b) Surface BZ for a cut nor¬ 
mal to y with projection of the “bent” BZ in red (thick) and 
a sketch of an occupied band pair along this loop. On the in¬ 
variant planes {k^ = 0 and = tt) the glide eigenvalue (color 
code) distinguishes the two bands. On the non-symmetric 
part (other values of k^) bands are generically nondegener¬ 
ate. The glide Chern-number is calculated by counting the 
winding number of </)(ks) following a single band (solid or 
dashed) around the X — T — X — M — Z — M — X loop, (c) 
Value of ^(ks) around the loop for the minimal model in the 
nontrivial phase. Color code indicates glide eigenvalue in the 
invariant planes. Note that following one band around, the 
phase winds 27r, indicating that the glide Chern-number is 
odd. 


the Berry-connection integral integrals in the y direction 

p2TT 

c^O^sit)) = / A^k)dky. (Bl) 

Jo 

where A^(k) = i (wkl Iwk) and the curve ks{t) = 
ikx{t),ky{t)) encloses half of the surface BZ (Fig[^(b)). 
In the case when there are multiple pairs of conduction 
bands with we naturally sum the Chern numbers for 
bands with the same glide eigenvalue. This quantity 
(mod 2) is a well defined topological invariant because 
pushing a Weyl node across one of the connecting sur¬ 
faces can only change the winding by 2. A definition that 
does not rely on distinguishing the two bands in the non- 
symmetric region is obtained realizing that the winding 
number (mod 2) is the same as the parity of crossings 
regardless of direction of an arbitrarily set branch cut by 
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throughout the invariant planes and one copy of the 
connecting plan^^. Even when the bands are degener¬ 
ate, phases of eigenvalues of the nonabelian Wilson loop 
operators along the y direction give equivalent quantities. 
This definition, while in principle well defined, still suffers 
from problems of distinguishing crossings from anticross¬ 
ings in a many-band model at finite k-space resolutioiP^. 

The fact that the Berry flux through a closed surface 
for a set of bands that are separated from all other bands 
on this surface is nonzero signals the presence a Weyl- 
node inside the bounded region connecting this set of 
bands to some other bands. This Weyl-node connects 
valence/conduction bands among themselves, the sys¬ 
tems we consider are fully gapped. This Weyl-node is 
also present in mirror symmetric insulators with differ¬ 
ent mirror Chern-numbers in the two invariant planes 


(sec. II A). We can think of the 3d band structure as an 
interpolation of 2d systems in the k^ky cuts as the tun¬ 
ing parameter kz evolves from 0 to tt. To interpolate 2d 
systems between the two cuts, we have to close and open 
a Dirac node between the occupied bands connected to 
different mirror sectors to transfer Berry flux, which cor¬ 
responds to a Weyl node inside the half 3d BZ. This Weyl 
node is topologically protected and must exist somewhere 
in the half BZ, mirror symmetry at the high symmetry 
planes prevents it from locally annihilating with its op¬ 
posite chirality mirror image that lives in the other half 
of the BZ. 


Appendix C: BdG formalism 

In this appendix we summarize general results of the 
BdG formalism used for superconductor^, with special 
emphasis to representations of symmetries in class C. 

The general form of a Hamiltonian for a superconduc¬ 
tor without any additional symmetry is 


At this point we digress to discuss transformation 
properties of the BdG Hamiltonian under physical sym¬ 
metry and gauge transformations. A general transfor¬ 
mation on the particle-hole degrees of freedom has the 
form 





Gk lEk 



(G3) 


because we require that the transformed particle and 
hole-like operators are still related by hermitian conju¬ 
gation. Moreover, preserving the fermionic commutation 
relations restricts the block matrix to be unitary. 

Now we turn to the case with spin rotation symmetry. 
We split the c vector in 2 halves for the spin z component 
Ck = ( Ck^ Ck^ ) and rewrite the Hamiltonian in a 4A^ x 
AN block form where N is the number of orbitals not 
counting spin. The U(l) spin rotation invariance around 
the z axis requires to commute with the generator of 
the rotations cr^r^ which restricts it to the block-diagonal 
form 


/ Ckt 

0 

0 

<5k \ 

0 

Ck; 

-<5^k 

0 

0 

-<51k 

-C^kt 

0 

V 4 

0 

0 



(G4) 


Imposing spin rotation symmetry for the x axis as well 
(commutation with a^Tz) means ^kt = Cki =: fk and 
(5k = If we now define a new set of operators dko- = 

^ Cko- ) with well defined momentum and spin 

quantum number, we can rewrite the Hamiltonian as the 
sum of two 2N x 2N terms for the two spin orientations 


H ^ dl.^H2kcrdkg 

kCT 


where 


(G5) 


■^ = 11(4 c-k 

k 




'fl2kt — 


Ck ^k \ 

4 ) 


(G6) 


= E (4 

where c is a vector formed of all the electron anihilation 
operators of the unit cell and the BdG Hamiltonian H\^ 
is a 2N X 2N block matrix with N orbitals in the unit 
cell. Ck = (Hermiticity) and Ak = —A^^. (Fermi 
statistics). It is customary to denote the Pauli matrices 
acting on the particle-hole space (the block structure of 
H) Ty. The paricle-hole symmetry is represented by C = 
TxJC for usual electronic systems such that = 1. This 
symmetry restricts the BdG Hamiltonian as 

Ci/kC-i = = -Elk (G2) 

which is automatically satisfied by this form of the Hamil¬ 
tonian. This transformation relates excitations with op¬ 
posite energy and momentum in the doubled spectrum. 


and i? 2 kt = ='■ I? 2 k with Pauli matrices 

acting on the space with the two components of d. The 
unitary relation between the two spin sectors guarantees 
that the spectrum is doubly degenerate, for every eigen¬ 
state there is another state with the same energy and 
momentum but opposite spin. 

The constraints on the form of i? 2 k can be summarized 
as 

Ci?2kC-i = ryH*_^,ry = -hl2k (G7) 

where we introduced the new particle-hole conjugation 
operator C = VyK, with = — 1 (in the body of the pa¬ 
per we use T instead of r for this set of Pauli matrices 
as well). This operator relates states with opposite en¬ 
ergy and momentum but the same spin. Note that this 
operator differs from the original particle-hole conjuga¬ 
tion in that it is combined with a spin flip, the physical 
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symmetry should also reverse spin. As the symmetry is 
antiunitary and squares to — 1 , a zero energy eigenstate 
at an invariant momentum must be doubly degenerate 
(on top of the spin degeneracy that is always present) by 
the same reasoning that proves Kramers degeneracy with 
= - 1 . 


Appendix D: Conventions for Bloch functions 


There are two widely used conventions to define the 
Bloch basis functions. When appropriate we use the 
convention where we define Bloch basis functions |x^) in 
terms of the orbitals of the unit cell |x^) = 10 ^) 

where R is the unit cell coordinate and a the orbital in¬ 
dex. Note the absence of phase factors corresponding to 
the position of the orbitals within the unit cell, so the ba¬ 
sis functions are strictly periodic in the BZ. While in this 
convention the information about the position of the or¬ 
bitals is lost, thus the polarizations computed via Berry 
vector potential integrals do not equal the true Wannier 
center positions, the Bloch Hamiltonian is BZ periodic, 
making some derivations more transparent. 


In the other convention we define |x^) = 

gjk(R-f.r„) 10 ^) ig position of the 

a-th orbital in the unit cell. The two conventions 
are related by the operator Wk with 
such that = Wk|Uk). Furthermore, the new 

Bloch wave functions obey the boundary condition 


l^k+c) = J and operators (including the Bloch 

Hamiltonian) satisfy Ok+G = where G is 

a primitive reciprocal lattice vector. This convention 
is usually assumed in formulae for electromagnetic 
response, as the naive Peierls substitution k —> k -f A 
only gives the correct phase factor for hopping in this 
case. The two conventions give equivalent results for 
quantized topological indices in most symmorphic cases, 
provided there is a continuous, symmetry preserving 
deformation of the lattice, such that all the orbitals 
are brought to the same point in the unit cell. In 
nonsymmorphic lattices however, this is never possible, 
as the shortest orbit of a point in the unit cell under 
the symmetry group modulo lattice vectors is longer 
than one, there is no crystal with one site per unit cell 
obeying a nonsymmorphic symmetry. For example with 
a single essential nonsymmorphic symmetry translating 
in the x direction one needs at least n lattice sites that 
can be arranged such that the positions are = txajn 
for a = l...n, so and Wg = I for 

perpendicular directions. In our glide examples this 
results in Wg^ = Pz, where p acts on the space of the 
two sublattices. 


Appendix E: Proof of quantization of 9 


In this appendix we provide a formal microscopic proof 
of our claim that any orientation-reversing space group 
(SG) operation quantizes 6. 

First we review the representations of space group op¬ 
erations in k-space. We use the convention with Bloch 
basis functions |xkO = X]r |‘/’r+x)> where we 

split the orbital index a = (x, Z), x = labels the sites 
of the unit cell by their real space position and I is an 
on-site orbital index accounting for spin, orbital angular 
momentum, etc. (the values I can take may depend on 
x). A useful property of this basis is that it is periodic in 
the real space coordinate (insensitive to the choice of the 

unit cell), i.e. / = |Xk*) for any lattice vector R. 


On the other hand, this basis is not BZ-periodic, instead 
k and k -f G are related by a constant transformation as 
explained in Appendix]^ We emphasize that our treat¬ 
ment is not specific to tight-binding models, the same can 
be told in the continuum, there x is the continuous index 
for position in the unit cell and I stands for the spin only. 
To go to the tight-binding approximation, we restrict the 
Hilbert-space to a finite set of orbitals per unit cell, the 
only assumption we make is that there exists a basis of 
localized states such that orbitals centered on different 
sites span orthogonal subspace^. 

Consider a general SG operation g = {0|t} (O is an 
orthogonal rotation and t a translation) acting on one of 
the basis states 



^0(R+x)+t 


(El) 


where U is the site and ^-dependent unitary represen¬ 
tation on the local orbitals, a double representation if 
the model is spinful. Applying this to the Bloch basis 
functions, with simple algebra we find 


5|Xk) 





(E2) 


with gx = Ox -|- t that is understood as a permuta¬ 
tion of sites at the same Wyckoff position. Grouping 
indices back together, this can be written as g |Xk) = 

^-^{Ok)tJJba 1^^^^ y(x,0,(x'.i') ^ 

The key observation is that in this basis the k- 
dependence decouples as a single phase factor. Con¬ 
sider the transformation of a Bloch eigenstate in the 
n-th band = ’^klXk)- symmetry transforma¬ 

tion results in a state at Ok, the coefficients transform 
as ( 5 n)ok = or in a compact notation 

9 (’^k) = ( 5 «)ok = As g is a symmetry op¬ 

eration, the transformed state is again an eigenstate of 
the Bloch Hamiltonian with the same energy, but at Ok. 
As a consequence, the transformation of occupied band 
projector operator Vk = J^ugocc. ''^kfik reads 


(57^)ok = ^ {gnhMl^ = ur^ul (E3) 

nGocc. 
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So if 5 is a symmetry, such that (<77^)k = "^k, any gauge 
invariant quantity that can be expressed through T^k are 
invariant if the k-space coordinates are also transformed. 

Now we are in a position to prove that 9 is invariant 
under orientation preserving SG operations and changes 
sign (thus becomes quantized) under orientation revers¬ 
ing ones. First consider a symmorphic SG operation 
g = {O|0} acting as (gn)ok = C/rik- Note that in this 
convention that reproduces the microscopic expressioiJ^ 
for the diagonal magnetoelectric coupling, the Berry 
connection is calculated using the coefficients only as 
^ram _ .^g drop terms that come from the 

derivatives of the basis states. It transforms under g as 


= i (5?T-)lc d(5TO)k = dC/wo-ik = 


(E4) 


— n d"™ 


The constant U cancels and the 3d BZ integral (§ for 9 
picks up a factor of det O from the point-group rotation 
of k-space, so g9 = (det O) 9. 

We will reduce the general case to the previous sym¬ 
morphic one. Consider a continuous family of trans¬ 
formations represented as (g (ki) n)ok = e“dOk t)fc 4 jj.^^ 
where t is fixed and 0 < A :4 < 1 is a tuning parame¬ 
ter, ^4 = 1 corresponding to the actual SG symmetry 
g{ki = 1) = g = {0|t}. We will show that 9 {k 4 ) cal¬ 
culated from a wave function transformed by g {k 4 ) is 
independent of k^. This is to be expected, considering 
the special case of 17 = 1 and 0 = 1 corresponds to a 
mere shift of the spatial origin by fc 4 t. We can calculate 
the change in 9 as ^4 changes from 0 to 1 by evaluating 
the second Chern form 1 ]^ with ki as the 

fourth coordinate. Using the expression (A4), as V is 


A: 4 -independent, we see that 9 is unchanged through this 
process, g (^4 = 0 ) has the same form as a symmorphic 
SG operation discussed above, those considerations are 
still valid even though g ( 0 ) is not a symmetry. 

The same result can be derived by direct substitution, 
but one encounters a subtlety we discuss now. Using the 
transformed wave function n'k = e**^*nk we find A'^ = 
substituting it in Q 


9' = 9-^ 

47r 




(E5) 


2tt. However, t can be an arbitrary vector if we interpret 
the above transformation as a shift of the entire crystal, 
or equivalently as a redefinition of the spatial origin. This 
shows that 9 is ill defined with nonzero Ghern numbers, 
we rationalize this observation below. 


This situation is analogous to the problem of polar¬ 
ization in 2d Ghern-insulator^^. From a mathematical 
point of view, there the polarization becomes ill defined 
because it is not possible to choose a gauge where the 
wave functions are BZ periodic in the presence of nonzero 
Ghern-number. Similarly, the Chern-Simons 3-form is 
not invariant under gauge transformations that are not 
periodic over the BZ, one needs to choose a periodic 
gauge to fix it s val ue, but the Ghern-number prevents 
this. Equation (A4) misses the extra correction term be¬ 
cause, while V stays BZ periodic throughout the defor¬ 
mation, A 4 does not, to apply Stokes’ theorem one has to 
add surface terms for 9BZx [0,1], that exactly reproduces 
the correction. However this does not make a difference 
in the cases with vanishing Chern-numbers discussed in 
this paper. Ai fails to be periodic because the bound¬ 
ary conditions for the Bloch functions change throughout 
the deformation. To guarantee that such surface terms 
do not appear during the deformation to the trivial state 
used to calculate 9, we prescribe that these deformations 
should be made keeping the lattice sites (thus the bound¬ 
ary condition) fixed, which is always possible. 


We can continue the analogy from a physical point 
of view. Polarization measures surface charge, but on 
the surface of a Ghern-insulator charge is no longer con¬ 
served due to the chiral anomaly. 9 measures the frac¬ 
tional part of the surface Hall conductance a^y. However, 
when the bulk has a Hall conductance per transverse unit 
cell (Jxyloiz^ depending on the definition of the “surface 
layer” it may contribute more or less to a^y. We may fix 
the boundary of the “surface layer” and push the crystal 
in the 2 : direction by G. The surface Hall conductance 
changes exactly by GCicy/a^, in agreement with our for¬ 
mal result. We conclude that there is no natural zero 
for 9 with nonzero Ghern-number, consequently our re¬ 
sult about the quantization of 9 is only meaningful if we 
restrict to the case of vanishing total Ghern-numbers, as 
we did throughout this work. 


To summarize, we proved that 9 transforms under a 
generic SG operation g = {0|t} as 


As t is constant the extra term can be expressed as a 
linear combination of Ghern numbers in various cuts of 
the Brillouin zone. If the Chern numbers vanish, the 
correction is zero, or if t is a lattice vector, it is an integer 
multiple of 27r that does not change the value of 9 modulo 


g9 = (det O) 9 (E 6 ) 

meaning that 9 is quantized to 0 or tt by any orientation- 
reversing SG symmetry. 




